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We combine the most recent observations of large-scale structure (2dF and SDSS galaxy surveys) 
and cosmic microwave anisotropics (WMAP and ACBAR) to put constraints on flat cosmological 
models where the number of massive neutrinos and of massless relativistic relics are both left arbi- 
trary. We discuss the impact of each dataset and of various priors on our bounds. For the standard 
case of three thermalized neutrinos, we find ^ < 1.0 (resp. 0.6) eV (at 2a), using only CMB and 
LSS data (resp. including priors from supernovae data and the HST Key Project), a bound that is 
quite insensitive to the splitting of the total mass between the three species. When the total number 
of neutrinos or relativistic relics A^eff is left free, the upper bound on ^ (at 2a, including all 
priors) ranges from 1.0 to 1.5 eV depending on the mass splitting. We provide an explanation of the 
parameter degeneracy that allows larger values of the masses when A'eff increases. Finally, we show 
that the limit on the total neutrino mass is not significantly modified in the presence of primordial 
gravitational waves, because current data provide a clear distinction between the corresponding 
effects. 



I. INTRODUCTION 

Neutrino properties are among the most difficult to be probed experimentally, due to the weakness of their inter- 
actions. Data from particle accelerators tell us that there are only three flavor neutrinos, while neutrino oscillation 
experiments show evidence for non-zero neutrino masses (for a recent review, see e.g. |lj]). Recent results strongly 
constrain the mass differences of the individual neutrino masses (actually masses squared, Am^) and mixing angles, 
but no definite conclusion can be drawn neither on the absolute scale of neutrino masses, nor on the existence of 
weakly coupled sterile neutrinos. Fortunately, cosmology is quite sensitive to the neutrino sector (see 01 for a review), 
and can shed light on these questions, as well as other interesting issues regarding the Early Universe, such as the 
process of neutrino decoupling from the primordial plasma. 

Currently, the most popular cosmological model is the flat adiabatic ACDM scenario, in which the present density of 
the Universe is shared between baryons. Cold Dark Matter (CDM) and a cosmological constant A. This model makes 
rather simplistic assumptions concerning the neutrino sector, consisting only of three ultra-relativistic neutrinos. It 
turns out that with a more refined description of the neutrino sector, one finds that only small corrections to the 
standard picture are allowed after comparing with current data on Cosmic Microwave Background (CMB) anisotropics 
and Large Scale Structure (LSS). However, these small corrections carry enough interesting physical implications to 
justify an active research effort, in particular after the first releases of the WMAP and SDSS data. The results 
of this effort are not only new cosmological bounds on neutrino properties but also a better understanding of how 
the errors depend on (i) the experimental CMB and LSS data, (ii) external priors on the cosmological parameters, 
(iii) intrinsic parameter degeneracies in the theory of cosmological perturbations, (iv) assumptions concerning the 
underlying cosmological model and parameter space. 

We here perform a new analysis using the most recent LSS (2dF, SDSS) and CMB (WMAP, ACBAR) data and 
an extended cosmological model with an arbitrary number of massive neutrinos and additional relativistic particles, 
parametrized via an effective number of neutrinos (Ncs). We extend the recent work of 3] and those that appeared 
after the release of WMAP data In particular, our underlying model is identical to that of ref. j^^li but our 

analysis differs since we include an extended set of data (such as the SDSS results and a more updated version of the 
2dF ones), a new prior on the matter density from SN-Ia |14j and non- linear corrections to the LSS power spectrum 
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on scales 0.1 h Mpc~ < k < 0.2 /i Mpc~ . Furthermore, we increase the number of free parameters to ten, in order 
to analyze the bounds in the presence of primordial tensor perturbations. 

The rest of the paper is organized as follows. After a short summary of the effects of neutrino masses and additional 
relativistic particles in Sec. ^1 we describe our analysis method and dataset in Sec. IIIII We discuss our results and 
compare with previous works in Sec. lIVI Finally, we conclude in Sec. IVl 

II. EFFECTS OF ADDITIONAL RELATIVISTIC PARTICLES AND MASSIVE NEUTRINOS 

Non-standard neutrinos and other weakly-interacting light particles leave their imprint on the evolution of the 
Universe, both at the level of background quantities and spatial perturbations. Here we describe the main effects of 
additional relativistic particles, massive neutrinos and their simultaneous presence. 



The density of radiation in the Universe is usually assumed to be given by that of photons and of three thermally 
decoupled neutrinos. These contributions are of the same order and fix the evolution of the Universe in the radiation- 
dominated epoch (RD). Thus, if the three neutrinos did not decouple thermally, or in the presence of sterile neutrinos, 
the total density of the Universe during RD (as a function of the photon temperature T^) would be significantly 
affected, producing a change in the time of equality between radiation and matter, and in the sound horizon at 
the time of decoupling. These changes are known to shift the angular scale of the acoustic peaks in the CMB 
anisotropy spectrum as well as their amplitude (mainly, through the early integrated Sachs- Wolfe effect). They also 
have an impact on the matter power spectrum P{k), because a shorter matter-dominated stage implies less growth 
for perturbations inside the Hubble radius. As a consequence, the wave-length corresponding to the maximum in 
P{k) will be shifted proportionally to the Hubble scale at the time of equality. Thus, the effect of N^g is mainly to 
change the background evolution. However, ultra-relativistic particles also have a smaller effect directly at the level 
of perturbations, explained in detail in Ref. [T5l |. 

All these effects can be parametrized by a single quantity: the effective number of relativistic degrees of freedom 
during RD, defined by the relation 



Here, Pr stands for the total energy density of radiation and = (tt^/IS) is the contribution of photons. The 
parameter A^eff is defined in such way that if neutrinos decoupled following the instantaneous decoupling approxi- 
mation, Nes just gives the number of flavor families. However, N^s could differ from three in the presence of extra 
relics (sterile neutrinos, light gravitons, gravitinos, majorons, effects from extra dimensions, etc.) or in the case of 
non-thermal decoupling. Actually, in the standard case a careful study of non-instantaneous neutrino decoupling 
shows that A^cff — 3.04 for three flavor families IT^. Note that A'cff is constant only when the neutrinos or the 
other relics are ultra-relativistic. 

The value of A'cff is constrained by Big Bang Nucleosynthesis (BBN) from the comparison with the measured 
primordial abundances of light elements. During the BBN epoch the nuclear reactions freeze out at a scale factor 
that depends on the expansion rate, which in turn is fixed by the total energy density of radiation. A BBN analysis 
shows that A'cff — 2.5^q'1) {2a) (see also J,8]), which is perfectly compatible with the number of flavor neutrinos. 

However, it is interesting to measure N^s independently of BBN (e.g. using CMB and LSS data) for at least two 
reasons. First of all, because the number Acff could change between the two epochs 0, |23|. A second reason is 
because the standard BBN model might be a good first-order description, but with possible corrections due to spatial 
inhomogeneities, leptonic asymmetries, etc., that could be evaluated with an independent measurement of A'cff. 



Neutrinos that possess masses larger or of the same order than the relevant photon temperature have different 
effects than a constant A'cff. For instance, neutrinos heavier than roughly 10~^ eV are not relativistic today. Neutrino 
masses have implications for the evolution of cosmological fluctuations, both at the level of background quantities and 
directly on the perturbations. 



A. Enhanced relativistic energy density (Ncs) 




(1) 



B. Massive neutrinos 
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It is well-known that massive neutrinos could account for a significant fraction of the total energy density of the 
Universe today, unlike relativistic thermal relics (fir ~ 5.6 x 10~^ per neutrino family). For fully non-relativistic 
flavor neutrinos, the contribution to the present energy density is directly proportional to the number density. For 
vanishing neutrino chemical potentials, the total neutrino contribution to the critical density is given by 



where h is the Hubble constant in units of 100 km s""'^ Mpc~^ and ^ runs over all neutrino mass states. For fixed 
neutrino masses, fli, would be enhanced if neutrinos decoupled with a significant dimensionless chemical potential 

= fJ-,y/T (or equivalently, for large relic neutrino asymmetries), simply because their number density would increase. 
In principle there exist some combinations of pairs (C(t'e)j C('^At,T)) that pass the BBN test and are not yet ruled out 
by the CMB+LSS limits on iVoff 21] . However, it was recently shown that the stringent BBN bounds on apply to 
all flavors, since flavor neutrino oscillations lead to flavor equilibrium before BBN 22, 23, 24] . The contribution of a 
potential relic neutrino asymmetry is limited to such low values that it can be safely ignored. 

When the density of the other fluids (photons, CDM, baryons, dark energy) is kept fixed, the sum over the neutrino 
masses ^ has a direct repercussion on the geometry of the Universe. If instead the spatial curvature is kept fixed, 
the total mass affects the relative contribution ilx of the other fluids. In any case, this background effect has an 
impact on the observable CMB and LSS power spectra. For masses of the order of 1 eV, this signature is rather small, 
but can be marginally detectable. 

In general, neutrinos tend to stream freely across gravitational potential wells, and to erase density perturbations. 
Free-streaming is efficient on a characteristic scale called the Jeans length, corresponding roughly to the distance on 
which neutrinos can travel in a Hubble time. For ultra-relativistic neutrinos, the Jeans length is by definition equal to 
the Hubble radius c/H, but for non-relativistic ones it grows at a slower rate than c/ H (in comoving coordinates, it 
even decreases with time during matter domination). Neutrinos with masses smaller than approximately Tdcc ^ 0.3 
eV are still relativistic at the time of last scattering, and their direct effect on the CMB perturbations is identical to 
that of massless neutrinos. For bigger masses, the decrease of the free-streaming scale is felt by perturbations which 
enter inside the Hubble radius before decoupling, which results in a small enhancement of the acoustic peaks with 
respect to the massless situation. 

In the intermediate mass range from 10^'^ eV to 0.3 eV, the transition to the non-relativistic regime takes place 
during structure formation, and the matter power spectrum will be directly affected in a mass-dependent way. Wave- 
lengths smaller than the current value of the neutrino Jeans length are suppressed by free-streaming. The largest 
observable wavelengths, which remain always larger than the neutrino Jeans length, are not affected. Finally, there 
is a range of intermediate wavelengths which become smaller than the neutrino Jeans length for some time, and then 
encompass it again: these scales smoothly interpolate between the two regimes. The net signature in the matter power 
spectrum is a damping of all wavelengths smaller than the Hubble scale at the time of the transition of neutrinos to 
a non-relativistic regime psi ] 



where n„i is the contribution of matter to the critical density. The damping is maximal for wavenumbers bigger than 
the current free-streaming wavenumber fcps 



We have summarized both the background and the direct effects of the neutrino masses on the CMB and LSS 
perturbations. The total signature is difficult to describe analytically. However, one should remember that for masses 
of order 1 eV or less, the dominant effect is the one induced by free-streaming on the matter power spectrum. 
Therefore, the usual strategy is to combine CMB and LSS measurements, where the former roughly fix most of the 
cosmological parameters, while the latter is sensitive to fcps and provides bounds on the neutrino mass. 




(3) 




(4) 



C. 



Combined effects 



In a situation with TV thermalized massive neutrino species, the cosmological model should include an equal number 
of parameters, namely (mi, . . . , itin)- However, at first order such a model could be described by only two parameters, 
N and the sum of all individual masses J2i=i n^^- This choice not only simplifies the problem, but also provides 
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the correct contribution of neutrinos to the total energy density both at early times, when all neutrinos are ultra- 
relativistic and the radiation density depends only on N, and at a late epoch when at least one neutrino mass is large 
compared with the temperature and the density is given in good approximation by the total mass. 

However, for a precise description TV and X]i=i n ™« only relevant parameters. As an example, let us 

compare two scenarios with three neutrinos but different mass spectra: a degenerate case with (mo, toq, TOq) and a 
case with (3mo,e,e), where e <ti ttiq. In both scenarios the neutrino density is the same for early and late stages 
of the Universe. But at intermediate temperatures of the order mo ^ ^ 3mo, the energy densities are different. 
It is easy to show that the expansion rate is temporarily enhanced in the second scenario, but this will only have a 
small signature in the CMB spectrum for mo ^ 0.2 eV. On the other hand, the free-streaming wavenumber fcps of the 
heaviest neutrino is larger in the second scenario. Thus, in principle we expect more damping and sharper bounds on 
the mass in the degenerate case. 

Nowadays, after many years of experimental effort, we know that neutrinos must be massive in order to explain the 
evidences for flavor oscillations from measurements of atmospheric and solar neutrinos, independently confirmed by 
data from the detection of neutrinos from artificial sources at experiments such as K2K and Kamland. These results 
lead to specific differences between the individual neutrino masses, at two different scales: Am^^j^ ~ 2.5 x 10""^ 
eV^ and Am^^^ ~ 7 x 10~^ eV^ But one of the masses remains unconstrained, which reflects the fact that 
oscillation experiments can not fix the absolute scale of neutrino masses. It is clear that, considering the present 
cosmological data, in order to have a measurable effect the three neutrinos should have roughly the same mass, 
following a degenerate scheme (mo, mo, mo), where mo ^ 0.2 eV. 

On the other hand, the positive results from the LSND experiment point to neutrino oscillations with Am^ 
0(1 eV^). The less disfavored scenario that could explain the LSND data, together with the results of atmospheric 
and solar neutrino experiments, contains 4 neutrinos following a 3+1 scheme where one of them is much heavier than 
the others (see for instance [2^) and the fourth neutrino must be sterile. It has been shown [23, that all four 
neutrino models lead to a full thermalization of the sterile neutrino before BBN, so they are disfavored by BBN. 
However, current CMB and LSS data can not completely rule out this possibility |7l[ia|. 

The most general scenario is that of a cosmological model with N^" thermally-decoupled massive neutrinos and 
extra relativistic degrees of freedom, parametrized by N'^ not necessarily integer. This model is described by a set of 
N'" + 1 parameters: (mi, . . . , m^fur, N"). However, since such a parameter space is too large for a systematic analysis 
we will consider restricted cases described by only two parameters: the total effective neutrino number iVeff = N'^ + N'^'^ 
and the total mass M. In order to check the impact of the distribution of the total neutrino mass among the individual 
states, we will study two cases: 

• the model that we call degenerate has A'cff massive neutrinos with the same mass. Let us emphasize that this 
is a simplified model where the physical interpretation of non-integer values of A'cft is not obvious. 

• a second model the we denote 1+r has only one neutrino with a significant mass, while the other N^s — 1 species 
are ultra-relativistic. 

These two models are not chosen arbitrarily. First, the degenerate model includes the standard situation with only 
3 flavor neutrinos degenerate in mass, while the 1+r model includes the 3-1-1 scenario described above. Second, and 
most importantly, the degenerate and 1+r models appear as limiting situations of the general case once the parameters 
A'off and M have been fixed. Indeed, the former has the smallest possible value of the free-streaming wavenumber, 
while the latter has the biggest /cps- For any intermediate model (like, for instance, the third case studied in [l3l |. 
with three effectively massless standard neutrinos and Ncs — 3 species with equal mass), the observational bounds 
deduced from the CMB and LSS observations should lay between those that we get in the degenerate and 1+r limits. 

III. METHOD AND DATA USED 

The WMAP spectrum and many other cosmological data can be accurately fitted with a six-parameter flat ACDM 
model described by the Hubble parameter h, the fractional density of matter Qm = 1 ~ ^A, the baryon density 

in dimensionless units lui, — flbh^, the optical depth to reionization r, and finally, the amplitude and the spectral 
tilt of primordial perturbations {Ag^ris). Most of our calculations correspond to a model with eight parameters: the 
six previous ones plus A^e£F and M, as previously defined for the degenerate and 1+r models. We will also study the 
consequences of the presence of a background of primordial gravitational waves, which would contribute to the CMB 
anisotropy spectrum. For this case, our parameter space will be ten-dimensional, adding the tensor-to-scalar ratio r 
and the tensor tilt nt. In all cases, we use the CMBFAST code 29| to calculate the power spectra. 

In order to compare the theoretical cosmological models with current observations, we use a Bayesian grid-based 
method described in some previous works (e.g. f^), instead of the widely used Monte Carlo Markov Chains (MCMC) 
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technique employed for instance in |12L |30j . The former has the inconvenient of being considerably slower from the 
computational point of view, and the advantage of being very robust if the hypersurfaces with the same likelihood 
in parameter space have very complicated shapes (for instance, an elongated and curved banana shape). Also, the 
MCMC method cannot deal with likelihood distributions with nearly degenerate local minima; mathematically, the 
existence of local minima cannot be disproved, but in practice, using current data and the wide variety of models 
discussed in the literature, such a situation never appeared. In order to analyze unusual models with many parameters, 
it is safe to use at least once a grid-based method, avoiding surprises that could arise from parameter degeneracies. 
However, our analysis always shows quasi-ellipsoidal likelihood contours for various combinations of parameters. Thus 
we believe that a MCMC method would give similar results. 

For each combination of data set and priors, we plot the Bayesian likelihood of each point in the two-dimensional 
space of neutrino parameters (M, A'eff) after marginalization over the other free parameters. Actually, for simplicity 
we approximate the integration process required for a true marginalization by a maximization of the likelihood over the 
remaining parameter space. Integration and maximization are known to be equivalent in the case of a multigaussian 
likelihood distribution and we expect the maximization technique to be reasonably accurate in our case, where we 
obtained quasi-ellipsoidal contours. We checked this explicitly by computing the marginalized likelihood of each of 
the six ACDM parameters: our results are very close to those of 0, il2j, obtained with a MCMC technique and 
marginalized by integration over the likelihood. Of course, our maximization routine needs to compute the likelihood 
not only at grid points, but also in between. At a given point in parameter space, this is done by interpolating 
quadratically each value of the C; and P{k) from the nearest neighbors of the grid. 

Our CMB data include the 1348 correlated points of the Wilkinson Microwave Anisotropy Probe (WMAP), which 
measure the temperature x temperature (TT) '311 ^^'^ temperature x E-polarization (TE) [3^ correlation functions 
on the CMB sky. WMAP provides the best available data on multipoles I < 900. For constraining the TT spectrum on 
smaller scales, we employ the results of the Arcminute Cosmology Bolometer Array Receiver (ACBAR) experiment [33, 
Is^ . We remove the highest band powers (probing the region / > 1800) which could be contaminated by foregrounds. 
The published ACBAR band powers are decorrelated, so it is self-consistent to use only the first 11 points of data. 

Our LSS data set consists of 32 correlated points from the 2 degree Field (2dF) galaxy redshift survey covering 
the range 0.02 h Mpc~^ < k < 0.15 h Mpc~^, and of 19 decorrelated points from the Sloan Digital Sky Survey (SDSS) 
[36| in the range 0.015 h Mpc~^ < k < 0.20 h Mpc~^. In order to compare the smallest wavelengths with the data, it 
is necessary to take into account small deviations with respect to the linear power spectrum. Following the analysis 
in [T^ , we compute non- linear corrections for each point in parameter space using the numerical procedure described 
in Appendix C of (33| . Since these corrections rely on N-body simulations of a pure CDM Universe they might not be 
perfectly optimized for hot plus cold dark matter. However, the models studied here include only a small fraction of 
non-relativistic neutrino density. Therefore, at first order the fitting formula of |3^ should be reasonably accurate in 
our case (and certainly better than introducing no corrections at all). Each redshift survey is expected to constrain 
the total matter power spectrum modulo a global normalization factor called the bias b. Unless otherwise stated, all 
our results were obtained after marginalizing over the 2dF and SDSS bias, treated like two free parameters. 

We will also use a prior on the current value of the Hubble parameter, measured by the Hubble Space Telescope 
(HST) Key project |33|: h = 0.72 ± 0.08 (Icr). Finally, we will impose some constraints on the current density of the 
cosmological constant deduced from the redshift dependence of type la supernovae luminosity. For a flat universe, 
Perlmutter et al. [s^ give a conservative bound Qm = 0.28 ±0.14 (Icr) that we denote as the SN99 prior. In addition, 
we will calculate the impact of using the more restrictive result = 0.28 ± 0.05 (Ict) from the recent work by Tonry 
et al. nil, that we label the SN03 prior. 

IV. RESULTS AND COMPARISON WITH PREVIOUS WORKS 
A. Degeneracies and priors 

In order to understand the impact of each data set and prior on the final results, we will introduce them step by 
step in the calculation of the likelihood. We will first focus on the degenerate model. 

We start using only CMB and LSS data without any prior, apart from the top-hat priors defined implicitly by the 
boundaries of our grid. We look for the two-dimensional probability distribution of the neutrino parameters in the 
range < Af < 2.25 eV and 1 < A'cff < 9. We find that most of this parameter space is allowed at 2a, as shown in 
panel (a) of Fig. ^ Actually some of the grid boundaries are reached by the allowed models , unlike in the rest of 
our analysis where they have no influence on the 2a allowed regions. 

The physical explanation for such loose constraints is well-known: when the density of relativistic relics is left free. 
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FIG. 1; Two-dimensional likelihood in {Ncs, M) space, marginalized over the six remaining parameters of the model. We plot 
the la (green / dark) and 2a (yellow / light) allowed regions. Here we used CMB (WMAP & ACBAR) and LSS (2dF & SDSS) 
data, adding different external priors as defined in section |llll (a) no priors, (b) HST, (c) HST-I-SN99, (d) HST+SNOS. 
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and iOm = ^mh^ ■ Indeed, one can vary these two quantities in the same proportion while keeping fixed 



'Cq; 



the 



redshift of equality between matter and radiation. In order to remove this degeneracy, it is necessary to impose some 
priors on h and Thus we repeat the same analysis, using now the HST prior. In this case, we obtain a band of 
aUowed models shown in panel (b) of Fig.^ This region stretches up to the maximal values of M and A'cfr in our grid, 
so we cannot derive yet limits on these parameters. However, the other grid bounds are now irrelevant, since none 
of the models allowed at 2a ever reaches them. The la preferred region includes models with (M = 0, A^cff = 3.04), 
showing that there is no evidence for massive neutrinos and/or extra relativistic degrees of freedom. However, it is 
interesting to see that large departures from the standard ACDM model cannot be excluded. 

These results show clearly how difficult is to put bounds on the neutrino mass in presence of an excess of relativistic 
relics during RD, and vice versa. For instance, if we assume A^eff = 3, we get a 2a upper bound M < 0.8 eV, while 
for TVgff = 6 this bound spectacularly increases to M < 2.2 eV. This trend was already observed in |3i IE 113 1 and 
can be explained as follows. Suppose that we start from the best-fit standard model with (M = 0, A^eff = 3), and 
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FIG. 2: Illustration of the main parameter degeneracy affecting our results. For three particular models, we plot on the left 
the CMB temperature spectrum (normalized to WMAP) and on the right the matter power spectrum (two of them have been 
rescaled by hand for clarity). See the main text for details. 

that we increase iVcfr in such way that the radiation density in the early Universe is multiphed by a factor a. Then, 
in order to keep the CMB power spectra roughly constant, we should maintain a fixed value both for Zcq and ^a. 
This would be achieved with a transformation of the type {ujrmh) — > {atOrm V^h). However, it is well-known that 
the shape of the matter power spectrum is given roughly by the parameter F = ^mh, that scales like \/aT with the 
previous transformation, meaning that the spectrum will have more power on small scales relatively to large scales. 
A neutrino mass can balance this increase through the free-streaming effect. Therefore, models with high values of 
A^eff are compatible with larger neutrino masses. 

We illustrate this parameter degeneracy in Fig. [3 where we show the CMB temperature spectrum (normalized to 
WMAP) and the matter power spectrum for three particular models. The first model (red/solid curves) is the best-fit 
model for the standard case {M,Nes) = (0,3). In the second model (green/dashed curves), we have increased the 
number of degrees of freedom to A^eff = 7 and performed the transformation on (wm, h) as previously discussed, leaving 
both Zeq and ^Im invariant. We have also increased ng a little bit. On the CMB figure it is difficult to distinguish 
the curves because they perfectly overlap, especially for the best measured scales (those of the first acoustic peak), 
while the difference at ^ < 10 is masked by cosmic variance. But as expected, the shape of the matter power 
spectrum is modified, with more power on small scales relatively to large scales. This leaves plenty of room for models 
with a significant neutrino mass and free-streaming effect, like the one featured here (blue/dotted curve) which has 
(Af,Aeff) = (2.25eV,7). 

In order to obtain bounds on the mass, it is necessary to add more restrictive priors. In particular, the supernova 
priors should help to better constrain flrn (around the central value 0.28) and therefore also the shape parameter F 
(around flmh = 0.28 x 0.72 = 0.20). This value of F is in very good agreement with the constraint obtained from 
LSS data, assuming no neutrino mass: the SDSS power spectrum points to F = 0.21 ±0.03 (2^. Therefore, we expect 
the SN prior to restrict the possibility of a large neutrino free-streaming effect and to improve the upper bound on 
M. This is what we observe by successively adding the SN99 and SN03 priors in our analysis. The corresponding 
contours, shown in panels (c) and (d) of Fig. ^ are the main results of this work. With the SN03 prior, our one- 
dimensional bounds on the effective neutrino number (marginalized over M) read 1.6 < N^g < 7.2 at 2a, in very good 
agreement with our previous results 1.4 < A'ofr < 6.8 5] obtained without SDSS, no SN priors and fixing M = 0. 
The one-dimensional upper bound on the total neutrino mass (marginalized over A'ofr) is M < 1.1 eV at 2a. For 
comparison, we list the bounds for fixed integer values of A'off are in Table n for different priors. Since we use more 
recent 2dF and SDSS data, as well as more restrictive SN priors, we are not surprised to find stronger bounds than 
Hannestad & Raffelt [l3|. 
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TABLE I: The 2a upper bound on the total neutrino mass M (eV), after marginahzation over the six cosmological parameters 
of the fiat ACDM model, for particular values of A'cff • We show the results for two limiting cases of splitting the total mass 
between the neutrino states: degenerate (all neutrinos with the same mass) and 1+r (one massive neutrino and the other treated 
as relativistic relics). Here we have used the full CMB and LSS data set, either alone (no priors) or combined with the HST 
and SN03 priors. For large A^cff values and in absence of priors, the upper bound is larger than the maximal value of M in our 
grid (2.25 eV). 

B. Role of LSS data 

So far we have used the LSS data as an indication of the shape of the matter power spectrum, but not its overall 
amplitude. This amplitude is difficult to measure, because of possible differences between the two-point correlation 
function of luminous galaxies and that of matter, a problem known as the bias uncertainty. The 2dF team has 
established that the bias b is almost scale-independent, and derived some constraints either on the redshift distortion 
parameter (3 = /h jU or directly on h These two results must be employed with great care since the bias 
is expected to depend on the mean luminosity and redshift of each particular galaxy sample. In order to use a 
self-consistent bias prior, we would need to compute some correction factors for each model (see for instance refs. 
[^|33,E3)- This technically difficult procedure, that relies on many assumptions, is beyond the goal of the present 
paper and we prefer to conservatively discard any bias prior, as in [H. Just for indication, we tried to repeat the 
previous analysis with a very naive bias prior. Instead of leaving the 2dF bias as a free parameter, we tried to add 
the constraint /b2d¥ = 0.43 ± 0.07 jlj to our full set of data and priors (which includes all the CMB-t-LSS data, 
the HST prior, and one of our two supernovae priors). As shown in panel (a) of Fig.|31 our results remain unchanged. 
This is consistent with the analysis of Elgar0y & Lahav who treat the bias prior in a detailed way and find no 
impact on the neutrino mass determination. 

Since the LSS data plays a crucial role in constraining the neutrino mass, it is worth comparing the impact of the 
2dF and SDSS power spectra. We go back to a data set consisting of CMB-I-LSS-I-HST-I-SN99, and remove either 
the SDSS or 2dF spectrum from the analysis. The results, shown in panels (b) and (c) of Fig.|3| should be compared 
with the combined analysis previously shown in (c) of Fig. ^ The SDSS power spectrum appears to be much more 
conservative, in good agreement with previous papers: the WMAP+SDSS constraint on the neutrino mass for A^off — 3 
is as large as M < 1.74 eV [13, while a WMAP-hother CMB-h2dF analysis gives M < 0.69 eV ||]. Consistently, our 
combined analysis gives intermediate results: for A^eff = 3 our WMAP-|-ACBAR+2dF+SDSS bound is M < 0.9 eV. 
Our results seem also consistent with the recent analysis in J^, where the corresponding bound M < 0.75 eV was 
found including data from both galaxy surveys up to A: ^ 0.15 h Mpc~^. 

Note that in order to employ the SDSS data until /cmax — 0.20 ft- Mpc~^, it is crucial to include the non- linear 
corrections to the matter power spectrum, in particular for the A^off bounds. The panel (d) of Fig. ^ was obtained 
with all the CMB and LSS data, plus the HST and SN99 priors, but in absence of non-linear corrections. A comparison 
with Fig.QJ; shows that the constraints on A'cfr are lifted by one unit. This can be easily understood. Fig. 0] shows a 
typical power spectrum with and without non-linear corrections. The linear power spectrum has less power on small 
scales, i.e. a smaller effective shape parameter F. As explained earlier in this section, this can be easily compensated 
by an appropriate increase in A^off, and h, while leaving the CMB spectrum almost invariant. 

C. Adding tensors 

We have discussed so far the dependence of our results under the choice of data set and priors. However, generally 
speaking, the bounds on a particular parameter also depend on assumptions concerning the underlying cosmological 
model, a simple six-parameter flat ACDM model in our case. It is clear that by adding extra physical ingredients 
that would compensate the effect of neutrinos, we would relax the bounds on (A^eff, M). It is inviable to perform a 
systematic study of all the ACDM variants proposed in the literature, but one of them deserves a particular interest. 
Indeed, the six-parameter ACDM relies on the existence of super-horizon cosmological fluctuations at early times, 
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M(eV) 

FIG. 3: Impact of various assumptions related to the LSS data set. The default model for which we show la (dashed) and 
2a (solid) contours was obtained with WMAP, ACBAR, 2dF, SDSS, plus the HST and SN99 priors. Each panel shows the la 
(green / dark) and 2a (yellow / light) allowed regions (a) when including the 2dF bias prior [41|. (b) without the SDSS data, 
(c) without the 2dF data, (d) without non-linear corrections. 

which strongly suggests that perturbations are of inflationary origin. But inflation also predicts a background of 
primordial tensor perturbations: the question is whether these gravitational waves are large enough to contribute to 
large-scale CMB anisotropics (fundamentally, this depends on the energy scale of inflation). It is thus important to 
see how the neutrino parameter bounds evolve in presence of two extra parameters, the relative amplitude and tilt of 
the primordial tensor s pec trum (r, rij). 

Previous analyses showed that for eight-parameter models (flat ACDM -I- tensors) a signiflcant contribution 

of gravitational waves is disfavored. However, one could expect that in a ten-dimensional model (flat ACDM -|- 
tensors, with Ncs and M), a new parameter degeneracy would show up and relax the various bounds. We performed 
such an analysis for our full CMB and LSS data set, adding the HST and SN03 priors. The resulting two-dimensional 
likelihood for (A'oft, Af ), marginalized over the other eight parameters, is shown in Fig.|Sl It is almost indistinguishable 
from that with a vanishing contribution of gravitational waves, which shows that the current data is clearly able to 
distinguish between the respective effects of tensors and neutrinos. Thus the cosmological bound on neutrino masses 
is robust with respect to tensors. This robustness also holds when including a non-adiabatic, incoherent contribution 
to the power spectrum such as those predicted by topological defects, as shown in a very recent work (Hj . 
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FIG. 4: Matter power spectrum for the best-fit model with (M, A^'cff ) — (0, 3) , plotted with and without non-linear corrections. 
The difference becomes significant in the region k > 0.15 h Mpc~^ probed by the last two points in the SDSS data set. 
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FIG. 5: la (green / dark) and 2a (yellow / light) allowed regions in (M, N^s) space, marginalized over the eight other 
free parameters of the flat ACDM-|-tensors model. For comparison, we show the contours corresponding to the case with no 
gravitational waves. The difference is very small, showing the ability of the data to make a clear difference between the effect 
of neutrinos and gravitational waves. 

D. Impact of mass splitting 

Our aim is to constrain cosmological models with an arbitrary number of massive neutrinos and with extra rela- 
tivistic degrees of freedom. However, so far we reduced the analysis to two parameters (M, Ngs), with the implicit 
assumption that all neutrinos were degenerate in mass. We have not discussed the fact that for a fixed total number 
of degrees of freedom N^s and total mass M, the evolution of cosmological perturbations depends on the splitting of 
the mass between the different species. 
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FIG. 6: la (green / dark) and 2a (yellow / light) allowed regions in (M, A^eff) space for the J -/-r model, compared with the 
degenerate model (dashed la and solid 2a contours). Here we have included all CMB and LSS data, plus the HST and SN03 
priors. 

As explained in sectional for fixed {M, N^s), in some sense the opposite case to the degenerate model is the 1+r 
scenario, for which all the mass corresponds to a single neutrino eigenstate, instead of being equally shared. We built 
a second grid of 1+r models and analyzed it with our most restrictive set of data (WMAP+ACBAR+2dF+SDSS) 
and priors (HST+SN03). We show in Fig.Elthe new allowed regions, compared with the previous ones. In the limit of 
small mass and small A^cff , the CMB and LSS power spectra are almost identical in the two cases, because the effect 
of the mass splitting is second-order with respect to that of the total mass. A priori, this does not guarantee that the 
(M, A'off ) iso-likelihood contours are asymptotically equal, because we are doing a Bayesian analysis and the contours 
are defined with respect to the best-fit model. For instance, in Ref. [l^ . the best-fit models for the degenerate and 
1+r cases correspond to high values of (M, Ncs) and are different from each other. This explains why the authors 
find different contours in the two cases even in the massless limit. The same occurs in our analysis when we do not 
impose any prior. However, when at least the HST prior is taken into account, the best-fit model is very close to 
(M, iVeff ) = (0, 3) and the contours only differ at high values of M and Nee- As expected from the physical discussion 
in Sec. ^ the model with only one massive neutrino is less constrained: remember that in that case, for a fixed M 
the free-streaming wavenumber fcps is larger. Therefore, the damping of the matter power spectrum is less efficient. 
The one-dimensional bounds on the effective number of neutrinos (marginalized over M ) now reads 1.6 < Ngg < 8.5 
at 2(7, and the limit on the total neutrino mass (marginalized over A'eff) increases to M < 1.5 eV {2a). The bounds 
for fixed integer values of A'off are given in the last column of Table U 



We have calculated cosmological bounds on neutrino masses and relativistic relics (A^cff) using the latest data on 
CMB (WMAP and ACBAR) and LSS (SDSS and 2dF galaxy surveys) in the framework of an extended flat ACDM. 
In the cases in which a comparison is possible, our results are in good agreement with those of previous analyses 
In the well motivated case of three flavor neutrino with degenerate masses, we found an upper limit on the 
total masses of M < 1.0 (resp. 0.6) eV using only CMB and LSS data and priors (resp. including priors on h and 
JIa). The bound for four thermalized neutrinos with only one of them carrying a significant mass is M < 0.8 — 1.2 
eV, depending on the priors used. Therefore, the 4-neutrino solution to the LSND results is not completely ruled out, 
but some tension with cosmological data exists, especially if the strong SN03 prior is taken into account. 

In the case of arbitrary A'off, our results are summarized in Fig. ^ and listed in Tabled They clearly show the 
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existence of a parameter degeneracy between the total neutrino mass and A^off, a trend already observed in previous 
works 0, B 0, 113 that we have explained in Sec. IIVI External priors on h and Q\ are found to be of particular 
importance for constraining respectively Ngg and M. 

Since the standard ACDM model (with its three effectively massless thermal relics) sits within the la preferred 
region, we find no evidence for exotic physics such as out-of-equilibrium neutrino decoupling, non-standard nucle- 
osynthesis, extra relativistic relics, a significant amount of hot dark matter, etc. However, large deviations from the 
standard case are still compatible with observations: for instance, a model with one neutrino of mass M = 1.5 eV 
and eight relativistic degrees of freedom is allowed by CMB and LSS data, even when all priors are included (see Fig. 
IS)). In order to exclude this model, it is necessary to take into account the prediction of standard BBN, which gives 
stronger limits on A'cff . 

The bounds obtained in this paper are based on the observation of cosmological perturbations (CMB and LSS), 
combined with constraints on the current expansion and acceleration rates of the Universe (HST and SN priors). 
Therefore they are completely independent from the predictions of primordial nucleosynthesis. It is remarkable that 
in the space of the two standard BBN free parameters (cj6,7Vcff), the preferred regions deduced from cosmological 
perturbations and from primordial abundances are perfectly compatible with each other, and more or less orthogonal: 
indeed, our analysis (with the most restrictive priors) gives 0.0215 < ojb < 0.0235 and 1.6 < N^s < 8.5, while standard 
BBN favors 0.017 < ujb < 0.026 and 1.6 < N^ff < 3.6 jl^] (all these bounds are at the 2a level). 

In order to test their robustness, we have also calculated the bounds on M and N^s in the presence of primordial 
tensor perturbations. Our results show that the bounds are practically unchanged, because current cosmological data 
is able to distinguish between the respective effects of tensors and neutrinos. 

Finally, we have considered the impact of a different splitting of the total neutrino mass among the individual states, 
an analysis also recently performed in ref. [T^ . We compared the case of complete mass degeneracy (all neutrinos 
with the same mass) with that where one neutrino state effectively possesses the whole mass. We found that the 
bounds on the degenerate case are more restrictive due to its more efficient free-streaming, in particular for larger 
values if Ncs. However, for three or four neutrinos the differences between the two cases are not significant. 

Our bounds are a clear indication that present cosmological data provide interesting bounds on the neutrino sector, 
complementary to those from terrestrial experiments. These include tritium beta decay experiments, which provide 
a current upper bound on the total neutrino mass of 6.6 eV at 2a while the KATRIN experiment 46] is planned 
to have an accuracy of the order 0.35 eV. Sub-eV sensitivity to neutrino masses is also expected for experiments 
measuring neutrinoless double beta decays (3^, but only for Majorana neutrinos and with a dependence on the 
details of the mixing matrix. 

However, the cosmological bounds should be taken with care, due to their dependence on the data (or priors) used, 
and also on the assumption of a particular underlying model. Examples are given by the works [3 El, where non-zero 
neutrino masses are preferred. This warning should not prevent us to be confident on the power of future cosmological 
experiments to limit (and eventually detect) neutrino masses and other neutrino properties. For instance, forecast 
analyses have shown that with future data there will be potential sensitivities to AA'off ~ 0.2 [TsL |49| (eventually 
improving BBN results) and neutrino masses of the order 0.1 - 0.2 eV with Planck and final SDSS data (HHHIsa, 
or with galaxy and CMB lensing ,5^ ^ . 
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